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Abstract 

The properties of baryonic matter have been investigated at finite density 
and temperature using different models. The variation of baryon masses and 
fractional number densities with baryon density and temperature obtained from 
different models have been compared. The quark hadron phase transition have 
been studied using Chiral Colour Dieletric (CCD) model in the quark sector. 
No phase transition has been seen for the different variants of the Zimanyi- 
Moszkowski model. However, a phase transition is observed for the linear and 
non-linear Walecka model. 







1 Introduction 



The study of hot and dense hadronic matter is a very old sphere of activity [|I|]. It 
is quite well known now that QCD is the fundamental theory of strong interaction 
where the fundamental particles are quarks and gluons. In principle, one should 
be able to study the low energy hadronic phenomena using QCD. But due to the 
nonperturbative nature of the problem, it becomes difficult to use QCD to describe 
hadronic matter. This leads people to calculate hadronic matter properties from 
effective models. 

The quark structure of hadrons lead people to believe that under extreme 
conditions, i.e. at very high temperature and/or density, the quarks, which are 
confined inside the hadrons, may have larger spatial extension 0, |^ compared to 
the typical size of a hadron which is about 1/m. So hadron to quark matter phase 
transition may occur inside neutron stars because of the very high density. On the 
other hand, a soup of quarks and gluons may be produced in heavy ion collisions. 
Such a matter was termed as Quark-Gluon Plasma (QGP) 0. Independent of the 
fact whether QGP will be produced, it can be said that a hot and dense matter 
will be produced under such extreme conditions. That can be either QGP or hot 
and dense hadronic matter. Furthermore, even if QGP is produced it will undergo a 
phase transition to hadronic matter. These phenomena make the study of hadronic 
matter an interesting of research. 

In this paper we have studied the properties of baryonic matter at high tem- 
perature and density from the five different hadronic models. These are the Walecka 
model (LW), the Zimanyi-Moszkowski (ZM) models and the non-linear Walecka (NW) 
model. The Walecka model |Q contains nucleons, scalar a-meson and vector w-meson 
and it reproduces the features of strong interaction i.e. short range repulsive and 
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long range attractive forces. Extensive studies have been carried out, after the in- 
vention of this model, towards the calculation of hadronic matter properties from 
this model 0. This revealed certain inconsistencies in the model. It was found 
that this model yields the incompressibility K = 524MeV compared to the accepted 
value K = 210 ± 30MeV. Furthermore the model gives a very low nucleon mass at 
the nuclear matter density. As a remedy of these problems other models of hadronic 
matter were proposed. In the NW model, in addition to the interactions present 
in the usual Walecka model, cubic and quartic self interactions of the cr-meson are 
also included. 

The ZM model differs from the usual Walecka model in the form of cou- 
pling between the scalar meson and the nucleon, which is a derivative coupling in the 
new model. As a result of this derivative coupling the model reproduces some of the 
experimental results nicely 0, §[. It yields the incompressibility K = 224A9MeV, 
which is much closer to the accepted value compared to the Walecka model, and a 
nucleon effective mass = 797.64MeV at the nuclear matter saturation density 
(p = 0.15/m~'^). But, the price one has to pay is that the model looses its renormal- 
isability due to the derivative coupling. This is acceptable as far as the discussion of 
the effective models goes, since the description of hadronic matter need only be valid 
upto the temperatures T < 200MeV provided the deconfinement phase transition to 
QGP is a reality. 

In the present work the LW and ZM models have been extended to include 
hyperons. Some efforts have already been made to include hyperon degrees of freedom 
in the nuclear equation of state ||, ^, [10|, [Tl|] but, most of these studies were restricted 
to zero temperature scenario. Section 2 is devoted to the study of different models 
of baryonic matter. In section 3 the qurak model, which is chosen to be the Chiral 
Colour Dielectric (CCD) model has been discussed. Section 4 contains the study 
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of the temperature and density dependence of baryon masses and number densities. 
The hadron-quark phase transition is discussed in section 5 and finally the results 
and conclusions are given in section 6. 



2 Hadronic Models 

We have five different hadronic models to describe the properties of baryonic mat- 
ter. These models are the Linear Walecka model, the three variants of the Zimanyi- 
Moszkowski model and the Non-linear Walecka model. The different variants of the 
ZM model is denoted by ZM, ZM2 and ZM3 (according to the notations followed in 
ref. |12|) models respectively. The ZM models discussed in refs. p], contain neu- 
trons and protons only and they have been used to describe symmetric nuclear matter. 
Here the ZM models as well as the linear Walecka model are extended to include hy- 
perons, p-meson and leptons (electrons and muons). The p-meson and leptons have 
been included to describe the asymmetric baryonic matter in /5-equilibrium. 

The Lagrangian density for the NW model is written as, 



-]-hmN{gaN(yf - jc{g„N(yY + H ^lih'^d^ - mi)ilJi (1) 

The Lagrangian densities fo LW model and ZM models can be written in a 
unified form as given below @, |^, [T2[| : 



i 
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^P;.Pr - l^yyj + T.M^7,^>^ - M,)^, (2) 



with a = 0, /5 = for the Walecka model, a = 0, /? = 1 for the ZM model, a = 1, /5 = 1 
for the ZM2 model, and a = 2, /5 = 1 for the ZM3 model. 

In the Lagrangians @) and @ mentioned above, ipi is the baryon field which 
can be neutron (n), proton (p), lambda (A) or sigma (S^) baryons, cj, p and a are 
the corresponding meson fields, g^i, Qui and Qpi are the couphng strengths of cr, and 
p mesons, respectively, with the corresponding baryons. The electrons and muons 
are denoted by ipi. The Non-linear Walecka Lagrangian includes cubic and quartic 
self-interactions of the a field with couplings h and c respectively 0. The quantity 
m* is given by 

-=(1 + ^) (3) 

where g^m is the coupling of scalar meson with nucleon, is the nucleon mass and 
(Jo is the mean field value of sigma field. 

In the Mean Field Approximation (MFA) the meson fields are replaced by 
their ground state expectation values. Then the equation of motion for the baryon 
field can be written as 

i^.d - M* - ^^,7ou;° + T37op3°] = (4) 

where T3 is the z-component of the isospin and M* is the effective mass of the baryon 
which is 

M; = M, - g.m*^ao (5) 
for the linear Walecka and ZM models. For the NW model M* can be written as, 

M; = M, - g^,a^ (6) 
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The mean field values of the meson fields are given by 



^0 = I]^(^i7oV'i) 



m,2 



Po = Ez:^(^^^'7o^«) (7) 



for all the models. The value of ctq for the Lagrangian is 

'm. 



mi I 
and for the non-linear Walecka model (|l]) is 

9< 



^^^^''^'^-^l (8) 



The Lagrangian has six parameters. These are gan/^a, guA/^a, gpn/^p 
and gp\/mp g^jn/fnuj and gi^A/rrii^. The nucleon couplings are determined from the 
nuclear matter saturation properties @, 12|. The values are given as follows: = 



357.4 and C^, = 273.8 for the Walecka model: 169.2 and 59.1 for the ZM model: 
219.3 and 100.5 for the ZM2 model and 443.3 and 305.5 for the ZM3 model where 
Cl = gl{MN/m^f and Cl = gl{MN/m^f. The value of c/p„ is 8.55451. The other 
three parameters are coupling constants of the hyperon-meson interactions and are 
not well known. These cannot be determined from nuclear matter properties since 
the nuclear matter does not contain hyperons. Furthermore, properties of hypernuclei 
do not fix these parameters in a unique way. Here the hyperon couplings are taken 



to be of the nuclear coupling ]13 



The non-linear Walecka model has eight parameters out of which five are de- 
termined by the saturation properties of nuclear matter. These are nucleon couplings 
to scalar (7o-/mo-, isovector gp/mp and vector mesons guijvti^ and the two coefficients 
in the scalar self interaction i.e. h and c. Here the p and lo meson masses are chosen 
to be their physical masses. 



The values of the parameters obtained are 



gan/m„ = 1.52505 x lO^^, g^^/m^ = 8.63433/782, gp/nip = 8.55451/770, 
b = 3.41799 X 10^3 and c = 1.46 x 10-^. 

Once again the hyperon couplings are taken to be | of the nuclear coupling. 



3 Quark Model 



The Colour Dielectric Model is based on the idea of Nielson and Patkos In this 



model, one generates the confinement of quarks and gluons dynamically through the 
interaction of these fields with a scalar field. In the present work, the chiral extension 
(CCD) of this model has been used to calculate the quark matter equation of state 
||13| . The Lagrangian density of CCD model is given by |]13| 



C = ^Pix){trdp - (mo + m/xix)U,) + il/2)g^^\aA';^ix)}^ 
+fjATr{d,Ud^U^) - l/2mj02(x) - {I / A)x\x){F;,{x)f 

+ {ll2)al{d,x{x)f - U{x) (10) 

where U = e*'*'"'^"/-^'' and f/5 = e^^'^'^"^^^^'' , ^{x), A^(x), x{x) and (/)(x) are quark, 
gluon, scalar ( colour dielectric ) and meson fields respectively, and m are the 
meson and quark masses, /^r is the pion decay constant, -F)jy(x) is the usual colour 
electromagnetic field tensor, g is the colour coupling constant and are the Gell- 
Mann matrices. The flavour symmetry breaking is incorporated in the Lagrangian 
through the quark mass term (mo + m/xU^), with mo = for u and d quarks. So the 
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masses of m, d and s quarks are m, m and mo + m respectively. The self interaction 
U{x) of the scalar field is assumed to be of the form 

U{x) = aBx\x)[l - 2(1 - 2/a)x{x) + (1 - 3/a)x\x)], (11) 

so that U{x) has an absolute minimum at x = and a secondary minimum at x = 1- 
Thus the parameter B is the analog of the bag pressure of the bag model. The 
interaction of the scalar field with quark and gluon fields is such that quarks and 
gluons can not exist in the region where x = 0. In the limit of vanishing meson 
masses, the Lagrangian of eqn. (|10]) is invariant under chiral transformations of quark 
and meson fields. 

The parameters of the CCD model are the quark masses rriq {q = u, d, s), strong 
coupling constant a^, the Bag pressure B and the constant a. These are obtained by 
fitting the baryon masses. In the present calculation we have used mg(^u,d) = 92MeV, 
TJiq^s) = 295MeV, 5^/^ = l52MeV, a = 36 and a, = 0.08. 

The quark matter calculation proceeds as follows. It is assumed that in the 
presence of non-zero quark/ ant i- quark densities, the square of the meson field ((</>^)) 
may develop a non-zero vacuum expectation value. This assumption is analogous 



to the assumption that, in the linear-sigma model ||T5[ the cr-field acquires non-zero 
vacuum expectation values. It has been found that this occurs when the quark density 
exceeds a certain critical value. As a result of this the effective quark masses decrease 
with the increase in (0^) and hence with increase in baryon density. Thus at large 
quark densities one obtains from CCD model an equation of state similar to the 
equation of state of free quarks and gluons. Furthermore, the colour dielectric field x 
and (0^) in quark matter are evaluated using mean field approximation. Quark-gluon 
interaction is incorporated upto two loop in thermodynamic potential. 
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4 Hot and Dense Baryonic Matter 

The thermodynamic equihbrium in the hot and dense matter can be imposed from the 
weak interaction conditions. This introduces the following relations among different 
bare chemical potentials as, 

f^A = l^n 
/^E- = l^n + IJ-e 
fJ'fi = fJ'e (12) 

In the medium the baryon chemical potentials are modified as, 

_ 1 

1 

fii:- = /^E- - gojs-i^o + 9pPo (13) 

where fti is the chemical potential of the i'th baryon in the medium and //j is the bare 
chemical potential for the same. 

The effective mass of a baryon in matter is as given in equations (11) and (12). 

The charge neutrahty gives, 

rip = rie + + n^- (14) 

In addition the baryon number conservation gives, 

riB ^ nn + np + nA + n-s- (15) 
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In the last two equations nt is the number density of the i'th particle and is 
given by 



r d^p 




1 (27r)3 





According to the usual convention the baryon chemical potential is defined as 

All the above equations need to be solved self consistently to obtain the den- 
sity dependence of baryon masses and number densities at different temperatures. In 
figures 1-3 the temperature dependence of n, A and S~ masses have been plotted for 
five different models which have been considered here and at three different densities 
- p = 0,0.15 and 0.3fm~^. The density dependence of baryon masses at three tem- 
peratures - T = 0, 100 and 200MeV have been plotted in figures 4-6. The number 
densities for different models at T = and 200MeV have been plotted in figures 7-11. 



5 Quark Hadron Phase Transition 

The quark-hadron phase transition can now be calculated by using the effective mod- 
els discussed in section 2 and section 3. The phase transition point is determined by 
adopting Gibb's criterion. This says that the point at which the free energies ( or 
pressure ) of the two phases, for a given chemical potential, are equal corresponds 
to the phase transition point. In the present context, the baryon chemical potential 
(Pb) is the only independent chemical potential. Then the crossing of the pressure 
curves for two phases in P -ps plane gives the phase transition point. Glendenning, 
on the other hand, considered a case where, in the mixed phase, neutron and quark 



matter are not charge neutral but the mixture as a whole is |16|. This aspect has 



been studied further by Heiselberg et.al. in ref [T^. This situation is not considered 
here. 
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For a first order transition, the derivatives of tlie P- fiB curve for tlie two 
phases at the phase transition point are not equal and the difference in the two 
derivatives gives the discontinuity in the density of the two phases at the transition 
point. The two phases coexist in this range of density. The latent heat of transition 
is given by the difTcrcncc in the energy densities of the two phases at the critical 
point. As mentioned earlier, the lattice QCD calculations indicate that hadron - 
quark phase transition may be weakly first order or second order. In a calculation, 
such as presented here, one would necessarily get a first order transition since two 
different models are employed to calculate the properties of quark and hadron phases. 

The pressure for Walecka and ZM models can be written as 

P = m*"^mlu', + m*"imj(pg)2 _ ^ ^ p^^^^^^ ^ ^ p^^^^^^ ^^^^ 

For the Non-Linear Walecka model pressure P is given by 

+ J2PFGirhi,fii) + J2PFGimi,IJ'i) (18) 

i I 

In the above two equations Ppc is the pressure for the free fermi gas and is 
given by 



For the quark sector the pressure is calculated upto two loop in the quark-gluon 
interaction and is given by 



(27r)'^ 

d^P . _g (r,)/T\ rriS-^ f 



16 2 [ (Ppnjp) „ f d^pd^q 



{2n)%{p) V {2n)%{p)e,{q) 

10 



X { 



2m 



*2 



+ 



(e<?(p) - eg(g))2 - e 



(e,(p) + e,(g))^ 

2 
3 



2 

- 1 



X 



n {p)n [q) + {p)n'^ (q) 
n~ {p)n~ {q) + {p)n~^ (q) } 



(20) 



In the above expression e^, eg and e,^ are the quark gluon and meson kinetic en- 
ergies, 7q, 7g and •j^ are the corresponding degeneracies, n~{p) = 1/(1 + exp{{eg{p) — 
jj)/T)) and n+(p) = 1/(1 + exp((eg(p) +/i)/T)) are the quark and antiquark distribu- 
tion functions respectively, N{p) = 1/ {exp{e^{jp) jT) — 1) is the distribution function 
for mesons and n{p) = n'^{p) + n^{p). In the quark sector the chemical equilibrium 
under weak decay and the charge neutrality give 



and 
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(21) 



-Uu - -Ud - -Us - He 



(22) 

The baryon density = | I]j '^i where i = u,d,s and the baryon chemical potential 
is defined as Hb = l^u + l^d + l^s- 

The P — fiB curves for different models have been plotted in figures 12-16. 



6 Results and Discussions 

In the present work we have studied the properties of hot and dense asymmetric 
baryonic matter in /^-equilibrium. The baryonic matter has been discussed within 
the framework of five effective hadronic models. They are the LW model [Q, three 
variants of the ZM model 0] - ZM, ZM2 and ZM3 and the NW model §]. We have 
extended the LW, ZM, ZM2 and ZM3 models to include hyperons and also discussed 
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the NW model. The models include four baryons - n, p, A and S~, the scalar a-meson, 
the vector mesons uj and p and the leptons (electron and muon) . The baryonic matter 
properties i. e. the masses and fractional number densities have been calculated from 
these models in the mean field level. The pressure of the baryonic matter has been 
calculated in the usual canonical fashion. The pressure for quark matter has been 



calculated from CCD model [13 



The temperature dependence of baryon masses at different densities have been 
plotted in figures 1-3 for different models. Fig.l is the temperature dependence of 
nucleon mass at different densities. This shows that the nucleon mass decreases with 
temperature and the maximum change in the nucleon mass is obtained in the LW 
model. Also, as the density is increased the baryon mass first increases slightly with 
temperature and then starts decreasing. The same qualitative behaviour is observed 
in the case of A and S~ masses (fig. 2 and 3). With density, the variation of nucleon 
mass is even faster. In the LW model the nucleon mass goes to zero at about O.Sfm^^ 
whereas in the other models the nucleon mass saturates in the range 400 — GOOMeV. 
The higher value is for ZM model and lower is for ZM3 model. 

In figures 7-11 the variation of fractional number densities, (nj/ns), of different 
particles with baryon density for different models have been plotted, each at two 
temperatures T = OMeV and T = 200MeV". Fig. 7 shows that the neutron number 
density in the LW model decreases and goes to zero at about 0.8fm~^. The proton 
number density increases initially from to 0.25/m~^ then it decreases slightly upto 
1.0/m~'^ after which it saturates. The A baryon appears at about ub = 0.25/m~^, 
increases till 1.0fm~^ to about 0.95 and then saturates. The S~ also appears at ub = 
0.25/m~'^, it increases first and then starts decreasing going to zero at ns = l.lfm^^. 
The matter becomes hyperonic at ub = 0.4/m~^. (We will use this terminology of 
"hyperonic matter", in the discussion of number densities, for the situation when the 



12 



A or S~ number density takes over n or p number density). In the same model, at 
T — 200MeV the matter is hyperonic right from the begining. The A number density 
varies from 0.39 to 0.85 as the density is varied from to 1.5/m~^. 

The nature of variation of number densities is somewhat different for the ZM 
model (fig. 8). At T = OMeV, the neutron number density decreases with baryon 
density but does not go to zero. The proton number density saturates at a higher 
value compared to that in the LW model. The transition from nuclear to hyperonic 
matter takes place at = l.l/m~^ which is considerably higher than that in the LW 
model. In the ZM model, at T = 200MeV, unlike the LW model, the matter is not 
hyperonic at the very begining. The transition takes place at about ub = 1.6/m~^. 
For ZM2 model (fig. 9), at T = OMeV, the transition takes place at ub — 1.6/m~^ 
whereas at T = 200MeV it is at = 1.9/m~^. Figure 10 gives an interesting result. 
This is for the ZM3 model. The matter does not become hyperonic ever. It stays as 
nuclear matter at all temperatures and densities. In the NW model at T = OMeV 
the transition occurs at = 0.85/m~^. Unlike LW model the neutron and sigma 
number densities do not go to zero. The proton number density also saturates at a 
higher value compared to that in the LW model. 

The hadron-quark phase transition has been studied using the Gibbs' criterion. 
The intersection of the two P — /ib curves, for two different models for the two 
sectors, gives the transition density. We have a first order phase transition due to 
the consideration of two different models in two phases. In figures 12-16 we have 
plotted the P — hb curves for five different hadronic models along with the quark 
model at five different temperatures. The quark-hadron phase transition, for this set 
of models, is observed only for LW and NW models and that too at T = OMeV 
and T = hOMeV. The transition, for the LW model at T = OMeV, takes place at 
liB = llOOMey which corresponds to Pc = 0.27/m~^. At T = bOMeV the critical 
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density is pc = 0.26/m ^. For the NW model these two values are pc = 0.47/m ^ 
and Pc — 0.37/m~^ respectively. 

To conclude, we have compared different hadronic models at finite temperature 
and density. These models give varying predictions for masses and number densities. 
This will have a strong bearing both for neutron star properties as well as heavy ion 
collisions in the laboratory. Also the neutrino emissivities will be different for these 
models. 

The hadron-quark phase transition has been investigated using CCD model. 
The phase transition is seen for LW and NW whereas none of the variants of ZM model 
gives the phase transition for the parameter set considered here, even at T = OMeV. 
This will imply the absence of quark-hadron phase transition inside neutron star. 
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